Analysis of transport properties of iron pnictides: spin-fluctuation scenario 
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We present a phenomenological theory of quasiparticle scattering and transport relaxation in the normal state 
of iron pnictides based on the simplified two-band model coupled via spin fluctuations. In analogy with anoma- 
lous properties of cuprates it is shown that a large and anomalous normal-state resistivity and thermopower can 
be interpreted as the consequence of strong coupling to spin fluctuations. The generalization to the supercon- 
ducting phase is also discussed. 
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I. INTRODUCTION 

Recently discovered superconducting (SC) iron pnictides 

|2l are at present in the focus of experimental and theoret- 
ical investigations in the solid state community. Besides high 
SC Tc the main motivation is the fact that a wide class of ma- 
terials opens yet another view on the interplay of magnetism 
and SC in metals. At present it seems that there is firm evi- 
dence that the physics of novel materials is not in the class of 
strongly correlated systems close to the Mott-Hubbard insula- 
tor as is the case in SC cuprates [3], investigated intensively in 
the last two decades. On the other hand there are challenging 
similarities, in particular in the anomalous transport proper- 
ties [0, m lal> in the presumable unconventional type of SC 
illlil, as well as in the importance of spin fluctuations appar- 
ently evidenced by the NMR relaxation |9] and the onset of 
spin-density-wave (SDW) order 1 10]. 

In this paper we put the emphasis on the transport prop- 
erties of iron pnictides (IP) and their understanding. Pub- 
lished experimental results on the d.c. electrical resistivity 
p{T) generally reveal high values at T > Tc, comparable in 
magnitude to underdoped cuprates At the same time, 

at elevated T ~ 300K the T-dependence is linear with large 
slope dp/dT, again similar to the well known anomalous vari- 
ation in cuprates. Since the systematics of different com- 
pounds is still not fully understood we concentrate here on 
the family emerging from the reference (undoped) compound 
LnFeAsO (LFAO) where besides the original Ln=La a vari- 
ety of other lantanides Ln=Ce - Dy has been investigated so 
far. The electron doping has been studied either by doping 
with F, i.e., in LnFeAsO i-a^F^; 10, El 01 or via oxygen defi- 
ciency LnFeAsOi_y fll2ll . The evidence so far is that different 
Ln do show similar results, while the resistivity p{T) reveals 
quite systematic and universal change with doping x [6] or 
y lfl2ll . E.g., the resistivity [13] changes from that of the 
SDW semimetal at x,y < 0.05 with large p{T > Tsdw), a 
property shared in particular by underdoped cuprates, over the 
intermediate (optimally doped) regime x,y ^ 0.1 with a quite 
pronounced linear law p (xT, into the overdoped regime with 
more Fermi-liquid (FL) behavior for y > 0.2. 

Also the thermopower S{T) is far from the FL behavior 



S cx T 10, Eh. reaching at T > Tc values characteristic for 
nondegenerate electrons, i.e., 15*1 ~ sq = fcs/eo = 86/iV/K 
again being a remarkable property of underdoped cuprates 
[13]. Similar message is emerging from strongly T-dependent 
Hall constant Rh{T) [14.] . So far, there are very few data on 
dynamical transport properties, nevertheless optical conduc- 
tivity in the same system seems to support non-Drude- 
like relaxation with large and cj-dependent transport relax- 
ation rate r^^ ^ uj. 

A detailed analysis of transport data on IP seems to be still 
premature due to mostly polycrystalline samples studied so 
far as well as due to the lack of doping systematics in the 
electron subsystem. Nevertheless in the following we argue 
that the similarity to cuprates arises from the strong coupling 
to spin fluctuations and the non-FL-like behavior following 
the marginal FL (MFL) scenario (l^. More specific origin 
of spin fluctuations and the spin-fermion coupling could be, 
however, quite different, e.g., due to the importance of the 
Hund's rule coupling Jh, also shown to lead to large p{T) 
inn ■ Also, emerging novel results on single-crystal transport 
obtained mostly by the electron doping the reference com- 
pound BaFe2As2 llisl 11911 require a quantitative reconsider- 
ation, not followed in this work in detail, due to observed 
higher mobilities. Nevertheless, several main characteristics 
as the linear p{T) law, large and anomalous thermopower 
S{T), T-dependent Hall constant Rh{T) again confirm our 
basic assumptions given below. 



n. A SIMPLIFIED MODEL 

A microscopic model for relevant electrons in IP appears to 
be quite complicated and still debated at present ]8]. While 
it has been argued that a (minimal) two-band model already 
contains the essential low energy physics in these materials 
1^0, 21, 22], the inadequacy of certain approximations to the 
effective low-energy band structure has been recently criti- 
cized ]23]. Nevertheless, for the purpose of our qualitative 
analysis we employ the model with two bands, one electron- 
like and another hole-like, coupled, however, through spin 
fluctuations introduced phenomenologically [24] and treated 
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within the lowest-order pertubation theory. It is evident that 
such a simpHfied model is not enough to describe calculated 
band structure 1I25I1 as well as the one observed via angle re- 
solved photoemission (26] or via the de Haas-van Alphen ef- 
fect |27], revealing (at least) four pockets at the Fermi sur- 
face in the actual Brillouin zone. However, in the simplest ap- 
proach to spin-fluctuation mechanism it is essential that SDW- 
type spin fluctuation inter-couple electron and hole bands, 
while the coupling among electron (hole) bands themselves 
should be less important. 

In contrast to IP, a vast experimental evidence in the last 
two decades shows that cuprates can be well modeled with a 
single-band model |3], e.g., the simple 2D Hubbard model or 
the t-J model [28] being, however, in the regime of strong 
correlations which can be only approximately described with 
a coupled spin-fermion model i24ll . Starting from the latter 
level assuming the MFL behavior of spin fluctuations and 
a strong coupling to electrons yield in cuprates anomalous 
quasiparticle (QP) damping ll29ll and transport relaxation llsoll 
as well as an unconventional SC | E9|] . 

We use a simplified model for IP describing the 2D system 
with one electron (eVband and the other hole (h) band crossing 
the Fermi surface ||8, 2^ 21]. I.e., in the (folded) Brillouin 
zone 10, E2I1 the h-like and e-like pockets are at fc ~ 0, and 
k ~ Q = (tt, tt), respectively. Within this effective model 
bands are coupled only via spin fluctuations, leading to 



kq,ss' 



,(1) 



and Ck, c?k^ (C^, C'') refer to electrons in e-like and h-like 
bands, respectively. We consider the corresponding Green's 
functions for e- and h-electrons 



where ^k — Ck ^ A^' cr = e,h{(T — h, e). 

III. QUASIPARTICLE DAMPING 



(2) 



Within the lowest order perturbation theory 1 2^ 3^ the self 
energies can be expressed as 

Sk(^) =32^771kq / / .912 ^ 



TT 

q - - 



cth 



(3) 



where Xq('^) is the dynamical spin susceptibility. 

To proceed we make several simplifications, which are 
expected to be a reasonable starting point for a qualitative 
analysis of transport quantities in IP. Spin response close 
to the antiferromagnetic (SDW) instability |10] centered at 



q '--^ Q = (tt, tt) we assume broad enough relative to h/e 
pockets to replace Xql"^) Xcii^) = x{'^)- In this case 
_ = E'^ and E(j - Ej; = E'\ and the QP damping is 
T'^iuj) = -ImE'"(tj). 

It should be however pointed out that the momentum de- 
pendence of the self energy could be important or even cru- 
cial. The latter would be the case if, e.g., the (accurate) nest- 
ing would play a role. Since anomalous behavior of IP seems 
to be quite robust, e.g., linear resistivity over a broad T and 
doping range, large thermopower S{T) as well as high Tc, we 
believe that this is not the case so that momentum dependence 
along the Fermi pockets is not crucial also requiring Xq(^) 
response not too narrow in q. 

With the above simplifications we get for the QP damping 

T%u) ^hj duj'giuo - J, uj')M%uJ - c^')x"(^'), (4) 

where rnq q — m-o,Q = m, A = w?, ■&x\A M'^ {yj) are the 
cr-band density of states (DOS). We also get 



^k(w) = -lm[0(w) - Ek + »r(w)]" 



(5) 



where = — ReE(a;) = ujZ^^{uj) defines the QP 

weight Z{uj). 

In analogy with cuprates, large and non-FL-like linear re- 
sistivity p{T) cx T and a presumable transport relaxation rate 
CX w in particular, require also a MFL behavior iHl of 
the spin fluctuation input x{^)- In om analysis MFL-type 
spin fluctuations are a phenomenological assumption which 
still has to be confirmed by specific experiments as, e.g., the 
inelastic neutron scattering. It is also fair to admit that a the- 
oretical understanding of such spin fluctuations in IP is not 
yet available, whereby even for cuprates there is not yet an 
agreement on the explanation for the MFL physics. Still, as 
for cuprates there are no accepted alternative scenarios for the 
non-FL behavior of QP damping and transport quantities. 

Hence, we employ the relation l28ll 



x"{io)=7rCiLu)thipLo/2), 



(6) 



where C{llj) — C{llj) + C{—lo) and C((jj) represents the 
dynamical spin correlation function. In analogy with the 
MFL scenario [1^ well established for cuprates we use a T- 
independent ansatz C{lo) Co, w < S> T |28, S^l- It 
should be noted that such an ansatz by construction satisfies 
the T-independent sum rule 13111 given by 







In underdoped cuprates the MFL behavior emerging from the 
above assumptions is qualitatively rather well established both 
experimentally |13l|32|] and from model calculations [28], al- 
though other forms close to Eq.® have also been proposed 
II32II . The origin of a non-FL behavior seems to emerge from 
the fact that due to the localized character of spins the low- 
u! spin-fluctuations exhaust the sum rule [31]. This is much 
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less evident for IP but still appears to be a prerequisite for 
the MFL-like behavior of transport quantities as described fur- 
theron. 

Assuming in the relevant regime J\f{uj) ~ A/" (being to low- 
est order unchanged even if the QP weight Z < 1) we get 

Viuj) = -f'^LUcoth{(3uj/2), = (37r/2)ACo7V'*, (8) 

i.e., approximately 

r'"(cj) ^ 7'"max(|w|,2T). (9) 

However, for further interpretation, in particular of the See- 
beck coefficient S{T), it seems essential that the DOS be 
nonsymmetric around w ~ 0. E.g., a possible assumption 
is J\f{uj ^ 0) ~ A/±, where the DOS can differ for w ^ 
and consequently 7'^ ^ 7 J . While the asymmetry of the re- 
laxation rate seems to be the only viable explanation for the 
large S{T), assumptions are expected to emerge from a more 
detailed analysis of Af{uj ^ 0) in a doped semimetal and cor- 
responding scattering rates for w ± 0. E.g., due to very shal- 
low hole bands at a; > for electron scattering in e-doped IP 
1 2^ 33t] as considered below one is effectively dealing with 
different A/±. 



IV. TRANSPORT QUANTITIES 

Turning to transport properties, we first consider optical 
conductivity cr(ci;) which we assume isotropic within the Fe- 
As, i.e., easy plane. Within the linear response (neglecting 
vertex corrections) a{uj) can be expressed as 13411 



k.cr 



where v^^ are the corresponding band velocities. In the fol- 
lowing we consider e-doped IP, therefore for simplicity we 
take as dominant the e-pocket contribution. Everywhere refer- 
ing to the e-band and defining the function 



N 



we get for slowly varying <I>(e) ~ $0 



cr(w) = $0 



1 - e-'^" 



dio 



,f{-Lo')f{u' -UJ)F, 



12 



"12 ^ ^12 



(11) 



(12) 



where F12 = r(w') + T{lo' - u) and = VL{u}') - - 

For a{uj ^ T) one gets from Eq. (fT2l i the extended Drude 
form 



a{iu > T) = $ 



(13) 



and 4> 



Sg/m*, with m* the QP mass in the 



e-pocket. Assuming further the MFL form for T{uo), Eq.®, 
the effective transport relaxation rate Ttri'jj) = 7qW where 
7a = Z(7++7_)/2. 



For the d.c. conductivity Eq.(fT2]i reduces to 1134,1351] 



1 



(14) 



which, on assuming constant $ and MFL-type F, Eq.®, im- 
mediately yields linear-in-T resistivity 



P 



T T 



An ^Aq ' 



(15) 



Within the same local approximation for Sk('jj) — '^{^) also 
the Seebeck coefficient S can be expressed as I35i] 



where 



S = —WSq, W = Ai/Aq, 



duj J 2Zr{uj) ' 



(16) 



(17) 



Under the MFL assumption for T{uj), all An are T- 
independent and in contrast to the FL behavior S oc T one 
gets a T-independent S ^ const. 

In the symmetric case 7_ = 7-(_ = 7 and by Eqs.([8]),(fT7]i 



Ao = 0.21/7, 



(18) 



while Ai — identically and therefore S{T) — 0. It is thus 
evident that a pronounced asymmetry in T{uj) is needed to 
explain large S{T) in IP, as discussed later. One situation 
possibly relevant for e-doped IP is that 7_ ^7+. This can 
happen if upon electron doping the h-pocket states diminish 
substantially at > as a source of scattering, leading to the 
reduction of 7+. In such a limiting case we get 



io = O.IO/7+, 



1.2 



(19) 



V. ANALYSIS OF EXPERIMENTAL RESULTS 



So far, most information for the normal(N)-state transport 
in IP is available for the d.c. resistivity p{T). We analyse 
here only the data for electron doped LFAO compounds. The 
whole range of x for F-doped LFAO has been measured re- 
cently li^, while S{T) as well as Hall coefficient Rh{T) have 
been measured also for x = 0.11 fl, H] and x = 0.05 ifTil . 
It should be reminded, however, that all data so far are for 
polycrystals while the theory is done for the transport oriented 
along the easy plane, so that measured p{T) should be at least 
scaled by some factor ^ < 1 to get the relevant planar resistiv- 
ity considered here. 

First, it is evident that inverse mobilities are larger in IP as 
compared to cuprates. Thus in LaSr2:CuOi_a; at T = 300K 
and for doping a; 0.03 - 0.1 |11], 



nhBoPab = (0.3 - 0.15)Vs/cm . 



(20) 
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Analogous results for IP compounds depend on the density 
X of carriers/formula unit, where rig — x/Vq and Vq is the 
volume of a formula unit. Even undoped x = OIP have finite 
but small n", i.e., xo > and it is plausible that x = xq + x > 
X. From data at T = 300K |6] we get 



fi ^ ~ neeop(T) ^ (5.6 — 8.1)a: Vs/cm 



(21) 



for X — 0.05 — 0.2, respectively. Similar results are obtained 
for X = 0.11 mi /Lt-i = 5.8iVs/cm2, and for x = 0.05 
1 1411 = S.OiVs/cm^. It is evident that such /i^^ are even 
higher than in low-doped cuprates [11]. 

Existing data for p{T) in LFAO show a large slope at higher 
T > Tc, also with a larger onset p{T ^ Tc) at x = 0.05, 
while for a: > 0.1 p{T < T* ^ lOOi^) becomes nonlinear 
and more FL-like. To estimate Aq and consequently 7 we use 
the slope at T < 300K, i.e., A^^ = ^dp{T)/dT and we get 



for a wide range of x | 

Ai) ^ av/x, 



a-3.8-10~^ 



7 ~ i,&x/v, (22) 



where v — m/m* . Similar values are obtained analysing 
other data for F-doped LFAO: for x = 0.11 iQ] and for 
X = 0.05 we get a 4.0 • 10"^ fl4]. Not much different 
is the development of p{T) for the oxygen deficient LFAO 
II12II where y < 0.03 compounds reveal a non-SC state with 
a SDW transition Tsdw > while for y > 0.1 again p{T) 
is nicely linear in T with a ^ 0.006. Assuming, e.g., val- 
ues from band structure calculations ~ 2 and x ^ 0.1 we 
arrive at very large 7 2.8. Such value 7 is most likely an 
overestimate due to too large p (all cited measurements are 
for polycrystals) and possibly due to relatively small QP mass 
enhancement v ^ 2. 

So far, there are only few data for optical conductivity 
(t(w). In the LFAO with As replaced by P (t((jj) was mea- 
sured and analysed f]3\ using the extended Drude fit yielding 
for Lj < lOOOcm^^ anomalous r^^ — Ttr ^ i.e., 7^1 
qualitatively consistent with the above estimates. 

Experimental results for S{T) for LFAO yield typically e- 
like 5 < with strong T-dependence with the maximum val- 
ues S ^ —So at T ~ lOOK. Assuming A7 = 7_ — 7+ > 7+ 
we get from Eq.lfTTIl 

w - A7/(127a) 

whereas in the extreme asymmetric limit, Eq.(fT9]l, we recover 
w = 1.2. 



VI. SUPERCONDUCTIVITY 

Let us finally comment on the relation of the above analysis 
to the treatment of the SC phase within the assumption of the 
spin-fluctuation induced pairing. We follow here closely the 
treatment of SC within the effective spin-fermion model 1291 
l36tl as derived from the microscopic strong-correlation model 
(planar t-J model) relevant for cuprates. Since our paper is 



not focused on the question of SC in IP, the aim is to connect 
parameters entering transport quantities to those determining 
the SC gaps and consequently Tc- 

For the discussion of SC equations can be generalized with 
Green's functions and selfenergies being 2x2 matrices Ii36ll 
whereby we (again) neglect e/h interband terms 

G£(c^) = Ko-e^Ti-S£H]-\ (23) 

In analogy to the normal state the lowest-order approximation 
for the self-energy can be written as 1 2^ 3^ 



S£(«t^n) = — ^TOkqG£_q(it^m)Xq(it^n"«W™), (24) 
q,m 

where iujn — «7r(2n + l)//3. Again neglecting the k depen- 
dence within each band, i.e., = Tl" one gets from Eq.Q 
a nonzero gap A'^ Z'^S^2(0)' 

-'^-,sc|5;£^.h^?s, ,.5) 



N ^ 
q 



where (E"^'^ 



2E?. 



h\2 



'q Q 



(A'*)^ and 



C| = I^{iuJn ~ 0)//q plays the role of the cutoff function 
with 



1 



(26) 



and/- =th(/3S-/2)/(2Sq 
Finally, at T = Eqs, 



reduce to 



A'^ = --Xx^Z^lS." 



de- 



(27) 



v/e2 + (A-)2' 

where u'^ are effective cutoffs. It is evident from Eq.dZTll that 
SC is of the s± - type [8], that is A'* = -r^A", > 0. As- 
suming also Af°'{e) ^ Af"' we get from Eqs.(l27ll.(l8]l. 



(28) 



with B = (4/7r) ln(ijJo/2r) connecting qualitatively the N- 
state transport parameters with the gap equation. Clearly the 
SC Eliashberg equations are treated in a simplified manner 
ll37ll in order to get the familiar BCS-type form. Still the 
message is quite clear: spin-fluctuation mediated interaction 
gives naturally the s± - type SC pairing consistent with other 
approaches [8]. Parameters entering Eq. (|28l l are besides 
B ^ 0{1) (depending on the form of the spin-fluctuation 
spectra) the cutoffs uj^ determined by the characteristic spin 
fluctuation frequency and 7°^. The latter have clearly the 
strongest influence and according to our estimates from trans- 
port, Eq.(l22b. 7°^ > 1 are large requiring the strong coupling 
approach both for the N and for the SC state. As discussed 
in the next section on the basis of emerging single-crystal re- 
sults, smaller resistivities p{T) are reported and consequently 
7 < 1. Still ^0 and the coupling remain at least moderate 
and in the same range as in optimum-doped cuprates, giving 
support for, or at least not contradicting, the notion for an SC 
pairing mechanism and its strength common to both IP and 
cuprates. 
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VII. DISCUSSION 



We have presented a theory based on the spin-fluctuation 
induced coupUng between the e- and h-bands in IP with the 
motivation to explain their anomalous N-state transport prop- 
erties. Existing experimental data on polycrystalline samples 
indicate that the QP damping and transport relaxation rates 
are even higher than in underdoped cuprates. It seems rather 
unlikely that quantitatively similar results should obtain for 
single crystals, as evidenced quite recently by measurements 
on, e^., BaFe2(Asi_xPx)2, a material from the 122 family of 
IP 113 811 or BaFe2-xCoxAs2 [19] where 7 and are substan- 
tially reduced with respect to values in (polycrystalline) LFAO 
samples of comparable doping. However, the marked linear- 
ity of p{T) ^ T/Aq + const, observed in BaFe2(Asi_xPx)2 
over most of the doping region with nonzero Tc testifies to 
the non-FL behavior, similar to cuprates. Likewise the evo- 
lution with doping in LFAO compounds is quite analogous, 
from a near insulator in an undoped substance to a FL-like 
behavior in the overdoped IP. Observe, however, that for the 
NdFeAsOi_xFx compound there is not much difference be- 
tween polycrystalline and single crystal data concerning p(T) 
and, e.g., a, Eq.(l22]i. Thus, for a; = 0.18 single crystal sample 
JH a rough estimate for T - 200K- 300K yields a - 0.023 
whereas in a polycrystalline sample x — 0.11 140] a ^ 0.009, 
resulting in 7 ^ 0.5 and 1.1, respectively, i.e., again compara- 
ble to cuprates. Thus a more systematic study of the transport 
and optical properties in single crystal compounds of the sep- 
arate families of oxypnictides is needed to settle this issue, 
particularly in view of the recent analysis of the competition 
between the linear and quadratic in T contributions to p{T), 
where the former is seen to dominate p(T) in sarnples with 
the highest Tc for the compounds there considered IsSH. 

Certain conclusions emerging from the above analysis still 
seem to be hard to avoid: a) the coupling to spin fluctuations 
is apparently substantial so that the QP damping is large with 
the QP overdamped in the low-doping regime, b) the effect 
of spin fluctuations on the N-state transport and on the SC 
pairing likewise appears to be strong implying pronounced 
spin fluctuations in the low-frequency window, both proper- 
ties shared by underdoped and optimally doped cuprates as 
well, c) if estimates emerging from experiments are correct 
the strength of the coupling could be beyond the applicabil- 
ity of the lowest-order perturbation theory employed here, d) 
the behavior of IP even at modest T > Tc is non-FL-like as 
clearly evidenced by large S{T) as well as the T-dependence 
of Rh{T), whereby the common features with cuprates stem 
from the large spin-fermion coupling and not from the Mott- 
Hubbard physics. However, rather scarce experimental evi- 
dence for low-energy spin fluctuations requires some caution 
and additional efforts to pin down the proper ingredients for a 
viable theory of IP. 



ACKNOWLEDGEMENTS 

Authors acknowledge fruitful discussions with C. Hess and 
the access to their unpublished data, and the financial sup- 
port of MHEST and JPSJ under the Slovenia- Japan Research 
Cooperative Program. T.T. acknowledges the support of the 
TRIP project. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. 

Chem. Soc. 130, 3296 (2008). 
[2] for an experimental overview see K. Ishida, Y. Nakai, and H. 

Hosono, J. Phys. Soc. Jap. 78, 062001 (2009). 
[3] M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 70, 

1039 (1998). 

[4] A. S. Sefatef a/., Phys. Rev. B77, 174503 (2008). 
[5] S. C. Lee et ai, J. Phys. Soc. Jpn . 78, 043703 (2009). 
[6] C. Hess etal., arXiv:081 LT60TV l. 

[7] K. Hashimoto et ai, Phys. Rev. Lett. 102, 017002 (2009). 

[8] for a theoretical overview see I. I. Mazin and J. Schmalian, 

arXiv:0901.4790 
[9] H. Mukuda et ai, arXiv:0904.4301 ^'1. 
[10] C. de la Cruz et ai. Nature 453, 899 (2008). 
[11] Y. Ando et al, Phys. Rev. Lett.87, 017001 (2001); Phys. Rev. 

Lett.93, 267001 (2004). 
[12] H. Eisaki et al., J. Phys. Soc. Jpn. 77, Suppl. C, 36 (2008); K. 

Miyazawa et ai, J. Phys. Soc. Jpn. 78, 034712 (2009). 
[13] J. R. Cooper and J. W. Loram, J. Phys. I (France), 6, 2237 
(1996). 

[14] Y. Kohama etal, Europhys. Lett.84, 37005 (2008). 

[15] M. M. Qazilbash etal., arXiv:0808.3748 i'l. 

[16] C. M. Varma etal., Phys. Rev. Lett.63, 1996 (1989). 

[17] K. Haule and H. Kotliai', New J. Phys. 11, 025021 (2009). 

[18] Y. Nakajima, T. Taen, and T. Tamegai, J. Phys. Soc. Jap. 78, 
023702 (2009)^ 

[19] L. Fang etal., arXiv:0903.2418V2; F RuUier-Albenque, D. Col- 
son, A. Forget, and H. Alloul, arXiv:0903.5243vl. 

[20] S. Raghu, X.-L. Qi, C.-X. Liu, D. J. Scalapino and S.-C. Zhang, 
Phys. Rev. B77, 220503(R) (2008). 

[21] W.-Q. Chen, K.-Y. Yang, Y. Zhou, and F-C. Zhang, Phys. Rev. 
Lett. 102, 047006 (2009). 

[22] M. M. Korshunov and I. Eremin, Europhys. Lett. 83, 67003 

(2008) . 

[23] S. Graser, T. A. Mayer, P. J. Hirschfeld, and D. J. Scalapino, 

New J. Phys. 11, 025016 (2009). 
[24] A. V. Chubukov and D. K. Morr, Phys. Rep. 288, 355 (1997). 
[25] D. J. Singh and M.-H. Du, Phys. Rev. Lett. 100, 237003 (2008). 
[26] C. Liu etal., arXiv:0806.2147v2. 
[27] A. I. Coldea et ai, Phys. Rev. Lett. 101, 216402 (2008). 
[28] J. Jaklic and R Prelovsek, Adv. Phys. 49, 1 (2000). 
[29] P Prelovsek and A. Ramsak, Phys. Rev. B65, 174529 (2002). 
[30] M. M. Zemljic and R Prelovsek, Phys. Rev. B72, 075108 

(2005). 

[31] P Prelovsek, I. Sega, and J. Bonca, Phys. Rev. Lett.92, 027002 
(2004); I. Sega and P Prelovsek, Phys. Rev. B 79, 140504(R) 

(2009) . 

[32] M. A. Kastner, R. J. Birgeneau, G. Shirane, and Y. Endoh, Rev. 

Mod. Phys. 70, 897 (1988). 
[33] Y. Sekiba et ai. New J. Phys. 11, 025020 (2009). 
[34] G. D. Mahan, Many-Particle Physics, Kluwer Academic 



6 



(2000). 

[35] G. Palsson and and G. Kotliar, Phys. Rev. Lett.80, 4775 (1998). 
[36] N. M. Plakida and V. S. Oudovenlio, Phys. Rev. B59, 11949 
(1999). 

[37] O. V. Dolgov, I. I. Mazin, D. Parker, and A. A. Golubov, Phys. 

Rev. B 79, 060502(R) (2009). 
[38] S. Kasahara et ai, arXiv:0905.4427V l 



[39] P Cheng, H. Yang, Y. Jia, X. Zhu, G. Mu, and H.-H. Wen, Phys. 

Rev. B78, 134508 (2008). 
[40] Z.-A. Ren et ai, Europhys. Lett. 82, 57002 (2008). 
[41] N. Doiron-Leyraud, P. Auban-Senzier, S. Rene de Cotret, A. 

Sedeki, C. Bourbonnais, D. Jerome, K. Bechgaard, and L. 

Taillefer, arXiv:0905.0964vl. 



